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ABSTRACT 



Context. The formation of planetesimals is often accredited to collisional sticking of dust grains. The exact process is unknown, as 
collisions between larger aggregates tend to lead to fragmentation or bouncing rather than sticking. Recent laboratory experiments 
have however made great progress in the understanding and mapping of the complex physics involved in dust collisions. 
Aims. We want to study the possibility of planetesimal formation using the results from the latest laboratory experiments, particularly 
by including the fragmentation with mass transfer effect, which might lead to growth even at high impact velocities. 
Methods. We present a new experimentally and physically motivated dust collision model capable of predicting the outcome of a 
collision between two particles of arbitrary masses and velocities. The new model includes a natural description of cratering and mass 
transfer, and provides a smooth transition from equal- to different-sized collisions. It is used together with a continuum dust-size 
evolution code which is both fast in terms of execution time and able to resolve the dust well at all sizes, allowing for all types of 
interactions to be studied without biases. 

Results. We find that for the general dust population, bouncing collisions prevent the growth above millimeter-sizes. However, if a 
small number of cm-sized particles are introduced, for example due to vertical mixing or radial drift, they can act as a catalyst and 
start to sweep up the smaller particles. At a distance of 3 AU, 100-meter-sized bodies are formed on a timescale of 1 Myr. 
Conclusions. Direct growth of planetesimals might be a possibility thanks to a combination of the existence of a bouncing barrier 
and the fragmentation with mass transfer effect. The bouncing barrier is here even beneficial, as it prevents the growth of too many 
large particles that would otherwise only fragment among each other, and creates a reservoir of small particles that can be swept up 

■ by larger bodies. However, for this process to work, a few seeds of cm in size or larger have to be introduced. 

Key words, accretion, accretion disks - protoplanetary disks - stars: pre-main-sequence, circumstellar matter - planets and satellites: 
O^l 1 formation 

f^*) [ 1. Introduction The micron-sized dust particles are coupled tightly to the sur- 

CN| ■ rounding gas, and their relative velocities are driven primarily 

■ One of the most popular planet formation scenarios is via core by Brownian motion. Since the resulting relative velocities are 
accretion, in which the formation of planets starts in the proto- sma ll, on the order of millimeters per second, the particles stick 

. ;_| . planetary disk with micron-sized dust particles that collide and together due to van der Waals forces. However, as the particles 

^ ■ stick together by surface forces, forming successively larger ag- increase in size, they become less coupled to the gas, and a num- 

H ! gregates (Mizuno 1980; Pollack et al. 1996). Traditionally, the ber of effects increase the relative velocities between them. For 

, _?3 . next stage in the formation process is the gravity-aided regime centimeter-sized particles, the predicted relative velocity is al- 

where planetesimals have formed that are so massive that the ready one meter per second, and meter-sized boulders collide 

gravity starts to affect the accretion and the strength of the body. at velocities of tens of meters per second. At these large col- 

However, to reach this regime, kilometer-sized bodies are lision energies, the particles tend to fragment rather than stick 

required, something which has proven difficult to produce due (Blum & Wurm 2008), which effectively prevents further growth 

to a number of effects such as fragmentation and bounc- (Dullemond & Dominik 2005; Brauer et al. 2008; Birnstiel et al. 

ing (Blum & Munch 1993), rapid inward migration (Whipple 2010). 

1972) and electrostatic repulsion (Okuzumi et al. 2011a,b). In the protoplanetary disk, gas pressure supports the gas 
A new planetesimal formation channel was introduced by against the radial component of the stellar gravity, causing it 
Johansen et al. (2007, 2011), in which mutual gravity plays to move at slightly sub-Keplerian velocities. Solid bodies do 
a role already between meter-sized boulders in turbulent and however not experience the supporting gas pressure, and instead 
locally overdense regions, resulting in a rapid formation of drift inward. As the particles grow larger, their relative veloc- 
kilometer-sized bodies. However, even the meter regime is diffi- ities compared to the gas increase, causing a significant head- 
cult to reach only by coagulation of dust aggregates. wind and a constant loss of angular momentum. At a distance of 
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1 AU, radial drift can cause meter-sized bodies to spiral inwards 
and get lost in the star on a timescale of a few hundred orbits 
(Weidenschilling 1977a; Nakagawaet al. 1986). 

These two obstacles give rise to the somewhat inaccurately 
named meter-size barrier (which ranges from millimeters to me- 
ters depending on the disk properties), above which larger bodies 
have difficulties getting formed. In order to reach the gravita- 
tional regime, bodies that are roughly nine orders of magnitude 
more massive are needed. 

The study of the dust evolution has however until recently 
primarily been done using simplified dust collision models in 
which colliding dust grains would either stick together or frag- 
ment (Braueretal. 2008; Birnstiel et al. 2010). The simplicity 
of the models has been a necessity because of large uncertain- 
ties and a small parameter space covered in mass, porosity and 
collision velocity in the laboratory experiments and numerical 
simulations. 

Recent years have however seen good progress in the lab- 
oratory experiments, as summarized by Blum & Wurm (2008). 
To provide a more complete and realistic collision model, 
Guttler et al. (2010) reviewed a total of 19 different experiments 
with aggregates of varying masses, porosities and collision ve- 
locities. In these experiments, the complex outcome was clas- 
sified into nine different types. Zsom et al. (2010) implemented 
this collision model in a Monte Carlo dust-size evolution code. 
The results showed clear differences from the previous collision 
models, and allowed for identification of the most important of 
the different collision types. They also found the important ef- 
fect of dust grain bouncing at millimeter sizes that halts the grain 
growth even before it reaches the fragmentation barrier. With the 
inclusion of a vertical structure, Zsom et al. (2011) still found 
the bouncing to be prominent, but also a number of other colli- 
sion effects caused by the vertical mixing. 

Progress has also been made with numerical simulations of 
dust (silica and ice) aggregate collisions using molecular dynam- 
ics codes (Wada et al. 2009, 201 1) with up to 10,000 monomers 
corresponding to aggregate sizes of around 100 yum. Based on 
these simulations, Okuzumi & Hirose (2011) developed a col- 
lision model where growth was possible for silicates up to ve- 
locities of 7 m/s, and for ices up to 70 m/s. By implementing 
relative velocities in the dead zone extracted from MHD simu- 
lations, they were able to form planetesimals made of ice, but 
not with silicates. Geretshauser et al. (2010, 2011) also devel- 
oped a dust collision code using SPH for particle sizes of cm 
and upwards. There is currently a discrepancy between the sim- 
ulations and the laboratory experiments, where the simulations 
have difficulties reproducing the bouncing events and generally 
observe much higher fragmentation threshold velocities. We will 
in this paper work primarily with the (more pessimistic) labora- 
tory data, but there is a great need to get the two fields to agree. 

One possible way to grow past the fragmentation barrier is 
so-called fragmentation with mass transfer, which was observed 
by Wurm et al. (2005) and can happen in a collision between a 
small projectile and a large target. The projectile is fragmented 
during the collision and a part of it is added as a dust cone to 
the surface of the larger particle, provided that the mass ratio 
between the two particles is large enough to avoid fragmenta- 
tion of the larger body. The mass transfer efficiency was studied 
by Kothe et al. (2010), who also showed that multiple impacts 
over the same area still lead to growth. Teiser & Wurm (2009a,b) 
have shown that growth of the target is possible even for collision 
velocities larger than 50 m s , and Teiser et al. (201 1) proved 
that the target could still gain mass even at large impact angles. 
These experiments have all shown that dust growth may proceed 



for large bodies at high velocities, and that this effect might even 
be able to produce planetesimals via collisional accretion. We 
discuss this process in more detail in Sect. 2.1. 

For the study of the dust-size evolution, the Monte Carlo 
approach of Ormel & Spaans (2008) and Zsom & Dullemond 
(2008) has a big advantage in that it allows for the simulation 
of a large number of particle properties and collision outcomes. 
A representative particle approach is used where a few particles 
correspond to larger swarms of particles with the same proper- 
ties. Each particle is given a set of properties, and each indi- 
vidual collision of the representative particles is followed. This 
approach uses very little computer memory, and adding extra 
properties cost very little in terms of execution time. If we want 
to study the effect of mass transfer, however, the Zsom et al. ap- 
proach has some problems, as it only tracks where the most mass 
is in the system. It therefore has difficulties resolving wide size 
distributions, which is required for the type of bimodal growth 
that the fragmentation with mass transfer effect would result in. 

Another method is the continuum approach, in which 
the dust population is described by a size distribution 
(Weidenschilling 1980; Nakagawa et al. 1981). The conven- 
tional continuum approach is the Smoluchowski method, where 
the interactions between all particles sizes are considered and 
updated simultaneously. This leads to very fast codes for a one- 
dimensional parameter-space (i.e. mass) compared to the Monte 
Carlo approach. Adding further properties such as porosity and 
charge is however very computationally expensive in terms of 
memory usage and execution time if one does not include tricks 
like the average-porosity scheme of Okuzumi et al. (2009). With 
the continuum approach, however, the dust is resolved well at 
all sizes, allowing for all types of dust interactions without any 
biases. It is also fast enough to follow the global dust evolution 
in the whole disk. 

The aim of this paper is to create a new collision model, de- 
scribing the outcome of collisions between dust aggregates of 
varying sizes and velocities, which is fast enough to be used with 
continuum codes. In this new model, we take into account the re- 
cent progress of the laboratory experiments, especially the mass 
transfer effect described above, and take a physical approach to 
transition regions from growth to erosion where the experiments 
are sparse. We then use this model in size evolution simulations 
using the local version of the code by Birnstiel et al. (2010) to 
study its implications for the formation of the first generation of 
planetesimals. 

The background of the new model and all the experimental 
work that it is based on is discussed in Sect. 2, and its implemen- 
tation is described in Sect. 3. In Sect. 4, we discuss the proper- 
ties of the disk in our local dust evolution simulations, as well 
as the implicit Smoluchowski solver that we have used. Finally, 
in Sects. 5 and 6, we discuss the results from the new model 
and show how the existence of a bouncing barrier may even be 
beneficial to the growth of planetesimals. 



2. Motivation for the new collision model 

Models to describe the growth of dust aggregates can generally 
be divided into two parts: A collision model describes the re- 
sult of a collision between two dust particles of arbitrary proper- 
ties (i.e. mass, porosity) and velocities. A dust evolution model 
uses the collision model to describe the evolution of the particle 
properties of an entire population of dust particles as they collide 
and interact with each other. In this section, we describe the lat- 
est laboratory experiments and our effort to produce a collision 



2 



F. Windmark et al.: Planetesimal formation by sweep-up 



model that can take these results into account while still stream- 
lining it to work well with continuum dust evolution codes. That 
means that the collision model cannot be as complex as the one 
developed by Guttler et al. (2010), but needs to focus on the most 
important collision types and aggregate properties. Still, we were 
able to include results which were not established or even known 
when the model of Guttler et al. was developed. 

2. 1 . Overview of recent experiments and simulations 

Numerous laboratory experiments have been performed to probe 
the collision parameter space of silicate dust grains, as summa- 
rized by Blum & Wurm (2008). This is a daunting task, as planet 
formation spans more than 40 orders of magnitude in mass and 
6 orders of magnitude in collisional velocity and collisional 
outcomes are affected by for example porosity, composition, 
structure and impact angle. The classical growth mechanism 
of dust grains is the hit-and-stick mechanism, which has been 
well studied both in laboratory experiments (Blum & Wurm 
2000, BW00) and in numerical simulations (Dominik & Tielens 
1997; Wada et al. 2009). Sticking collisions are also possible via 
plastic deformation at the contact zone (Weidling et al. 2012, 
WGB12) and geometrically by penetration (Langkowski et al. 
2008, LTB08). 

Previous collision models have due to limited data 
with few exceptions only included sticking, cratering and 
fragmentation with simplistic thresholds (Nakagawa et al. 
1986; Weidenschilling 1997; Dullemond & Dominik 2005; 
Tanaka et al. 2005; Brauer et al. 2008). In order to study the 
effect of the progress in laboratory experiments, Guttler et al. 

(2010) and Zsometal. (2010) presented a collision model 
containing nine different collisional outcomes and used this 
in the Monte Carlo dust evolution code developed by 
Zsom & Dullemond (2008). Their model contained three addi- 
tional types of sticking collisions besides the normal hit-and- 
stick, and they also identified two growth-neutral bouncing ef- 
fects and three different fragmentation effects in which the 
largest particle is eroded. It was found that several of the new 
collision types played a role for the dust-size evolution, which 
proved the necessity for a more complex dust collision approach 
than what had been previously used. Before even reaching the 
fragmentation barrier, at which fragmentation events between 
similar-sized particles prevent further growth, they found the ex- 
istence of the so-called bouncing barrier. Bouncing collisions be- 
tween smaller particles of intermediate velocities proved to be an 
efficient barrier for growth already at millimeter-sizes. It should 
be clarified that bouncing is, in principle, not bad for growth. 
What makes the bouncing barrier a problem is the lack of stick- 
ing over such a large range of masses and velocities that there is 
no way for the particles to grow further. 

Bouncing between dust aggregates is at the moment a 
hotly discussed topic. It has been reported from a large num- 
ber of laboratory experiments of different setups and mate- 
rial properties (Blum & Munch 1993; HeiBelmann et al. 2007; 
Langkowski et al. 2008; Kelling & Wurm 2009; Guttler etal. 
2010; Weidling et al. 2012), but molecular dynamics simulations 
show significantly less or no bouncing (Wada et al. 2007, 2008, 
2009; Paszun & Dominik 2009). These rebounding events hap- 
pen in collisions where the impact energy is so high that not all 
can be dissipated by restructuring of the aggregates. Wada et al. 

(201 1) argue that this would happen only for very compact ag- 
gregates where the coordination number is high, which is in con- 
tradiction to what is seen in the laboratory. We will in this work 
base our model on the laboratory experiments, but this is a very 
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Fig. 1: The size-size parameter space of the dust collision exper- 
iments (blue boxes) laying basis for the new collision model. All 
marked experiments are discussed in more detail in Sect. 2.1. 
The contours and gray labels mark the collision velocities in 
cm s~' expected from the disk model described in Sect. 4.1, and 
do not always coincide with the velocities studied in the experi- 
ments. 



important matter for the dust growth and will need to be investi- 
gated further. 

In our new model, we implement the most important colli- 
sion types identified by Zsom et al. (2010), and also take into ac- 
count the results from a number of recent experimental studies. 
Many new experiments have been performed that have further 
increased the understanding of the collision physics of dust ag- 
gregates. In Fig. 1, we plot the parameter space of a selection of 
important laboratory experiments laying basis for the new colli- 
sion model. 

Provided that the mass-ratio between the two involved parti- 
cles is large enough (from now on called the projectile and the 
target for the smallest and largest particle, respectively), the pro- 
jectile can fragment and parts of it stick by van der Waal forces 
to the surface of the target. This was studied by Wurm et al. 
(2005) and Teiser & Wurm (2009b, TW09b) for millimeter to 
centimeter-sized projectiles shot on a mounted decimeter-sized 
dust target at velocities up to 56.5 m s _1 . It was found that the 
accretion efficiency even increased with velocity, and could be 
as high as 50% of the mass of the fragmented projectile, where 
Guttler et al. (2010) only assumed a constant 2%. The effect was 
also observed by Paraskov et al. (2007, PWK07) in drop tower 
experiments where also the target was free-floating without a 
supported back. The mass transfer efficiency at slightly lower ve- 
locities (1.5 - 6 m s _1 ) and with millimeter-sized projectiles was 
studied in more detail by Kothe et al. (2010, KGB 10), who con- 
firmed the velocity-positive trend. Teiser & Wurm (2009a) and 
Kothe et al. also studied multiple impacts over the same area, 
and could conclude that growth was possible even then, without 
the newly accreted material being eroded. It was also found that 
growth was possible even at very steep impact angles. Beitz et al. 
(2011, B+ll) performed experiments between cm-sized parti- 
cles at even lower velocities (8 mm s _1 to 2 m s~') and found 
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mass transfer events already at 20 cm s _1 , right at the onset of 
fragmentation. 

This mode of growth, where small projectiles are impacting 
large targets, is related to the work of Sekiya & Takeda (2003, 
2005). They did numerical studies on whether small fragments 
formed in an erosive collision could be reaccreted back onto the 
target due to gas drag. The conclusion was that if the fragments 
were /mi-sized, and the target large enough, the gas flow would 
actually lead the fragments around the target and thereby pre- 
vent reaccretion. For the mass transfer effect, it is important to 
verify whether this effect could prevent 100-1000 yum-sized pro- 
jectiles from impacting on the target in the first place. The im- 
portance of this effect can be estimated with a simple comparison 
of timescales using reasonable parameters from the disk model 
(discussed in more detail in Sect. 4). The stopping time of a small 
particle is given by 

t s = ~ 2500 s (1) 
Pg ■ 11 

where £ ~ 1 g cirT 3 is the solid density of the projectile, a p ~ 
100 //m is its radius, p g ~ 10~ 10 g cm -3 the expected midplane 
gas density at 3 AU and u ~ 4 ■ 10 4 cm s the mean thermal 
velocity of the gas. The time it would take for the projectile to 
pass the target is given by 

T"pass = ~ °- 2 s (2) 

where a t ~ 100 cm is a typical target size and Av ~ 5000 cm 
s _1 the relative velocity between the particles. Since t s >> T pas s, 
it would take too long for the projectile to adjust to the gas flow 
around the target, and the two particles would collide. If the pro- 
jectile was instead only 1 p,m in size, the timescales would not 
differ so much, and the gas flow might play a role. 

Another recent experimental progress is the refinement of 
the threshold velocity for destructive fragmentation, where the 
target is completely disrupted. Beitz et al. (201 1) performed ex- 
periments to determine the onset of global fragmentation of the 
particles, and found that cm-sized particles fragmented already 
at 20 cm s , much below the lms" 1 threshold found for mm- 
sized particles by Blum & Munch (1993). This points towards a 
material strength that decreases with mass, as predicted for rocky 
materials among others by Benz & Asphaug (1999). This can 
be explained by a larger probability for faults and cracks in the 
material the larger the particle, and it is along these cracks that 
global breaking and fragmentation takes place. No experiments 
have as of yet been performed to study the fragmentation thresh- 
old between different-sized dust aggregates, but one can gener- 
ally assume that the velocities needed would be higher with an 
increasing size-ratio, as is seen both in experiments and simula- 
tions of collisions between rocky materials (Stewart & Leinhardt 
201 1; Leinhardt & Stewart 2012). 

In order to provide more data in the transition region be- 
tween sticking and bouncing, Weidling et al. (2012) studied col- 
lisions between particles of 0.5 - 2 mm in size and at velocities 
of 0.1 - 100 cm s _1 . In these experiments, sticking collisions 
were found (in coexistence with bouncing events) for higher 
velocities than was previously expected (Blum & Wurm 2000; 
Guttler et al. 2010), and enough data now exists to define a tran- 
sition regime between only sticking and bouncing. Similar ex- 
periments with smaller particles roughly 100 yum in size were 
performed by Kothe et al. (unpublished), and are consistent with 
the threshold of Weidling et al. 



Schrapler & Blum (2011, SB11) have also performed ero- 
sion experiments between /mi-sized monomer projectiles and 
mounted high-porosity aggregates for velocities up to 60 m s _1 , 
in order to determine the erosion efficiency as a function of the 
collision velocity and the surface structure. They found the ini- 
tial stages of the monomer bombardment to be very efficient 
even at low velocities, but after the most loosely bound monomer 
chains had been knocked off and the surface had compacted, the 
erosion was greatly decreased. 

2.2. Individual treatment of collisions 

In the collision model of Guttler et al. (2010), a binary approach 
was used for the particle mass ratios and porosities. Below a 
certain set critical mass ratio, r c = m t /m p = 10, 100, 1000, 
the collision was treated as being between equal-sized particles, 
leading for example to global fragmentation if two large parti- 
cles collided at high velocities. If the mass ratio was above the 
critical ratio, the particles were treated as different-sized, and a 
high- velocity collision would instead lead only to cratering. The 
same approach was taken for the porosity. Below a critical poros- 
ity (p c = 0.4, a particle was considered porous for the purpose 
of determining the collision outcome, and above it, the particle 
was considered compact. Combining these two binary proper- 
ties gave eight different collision scenarios, where the collision 
outcome was determined by the particles masses, porosities and 
relative velocity. 

In the new model, we instead use the current laboratory data 
to do an interpolation between the two extreme mass-ratios. This 
provides a continuous transition from equal-sized to different- 
sized collisions, and allows us to distinguish between collisions 
of different sizes at intermediate mass-ratios, and provides a nat- 
ural and smooth transition between the two extremes. We can 
therefore obtain what velocity is needed to cause global frag- 
mentation for a specific mass-ratio, which gives us a more pre- 
cise tool to study the effects of where global fragmentation turns 
into local cratering. 

It is however necessary for us to make a simplification re- 
garding the porosity of the dust grains. Adding additional prop- 
erties to the dust grains is very computationally expensive for 
continuum codes like the Smoluchowski solver that we use for 
the dust-size evolution, compared to Monte Carlo codes. In the 
Monte Carlo approach, each timestep only consists of one col- 
lision between a representative particle and a swarm of identi- 
cal particles. After the collision, the properties (i.e. mass, poros- 
ity, charge) of the representative particle is updated, and a new 
timestep is initiated. This means that for a simulation with n rep- 
resentative particles, each new property only adds an additional 
time 0{n) to the execution time. 

In the Smoluchowski method, one has to numerically solve a 
number of differential equations to update the number density of 
all mass bins. For each grain size, n 2 interaction terms need to be 
considered, where n is the number of mass bins. This is because 
a mass bin can collide with all bins including itself, but frag- 
mentation can cause mass to be put into it also from a collision 
between two other bins. If an additional property such as poros- 
ity is included, m — n porosity bins would need to be included. 
For each n ■ m bin, (n ■ m) 2 interactions would now need to be 
considered, and the code would be slower by a factor 0(m i ). In 
order to include porosity in the Smoluchowski solver, we would 
therefore require some analytical trick like an average porosity 
for each mass bin described in Okuzumi et al. (2009). This is 
however outside the scope of this paper, and we instead assume 
that all particles are compact at all times. This is likely a good ap- 



4 



F. Windmark et al.: Planetesimal formation by sweep-up 



proximation for larger particles outside the hit-and-stick region, 
as bouncing collisions quickly lead to compaction of the parti- 
cles. This finally gives us one single collision scenario, where we 
can for a collision between any two given particles determine the 
outcome based on their masses and relative velocity. 



3. Implementation of the model 

We discuss in this section the details on how the new colli- 
sion model has been created and implemented into the code. 
We choose to include only the collision types that proved to 
be most important in the simulations of Zsom et al. (2010). The 
collision types considered here are sticking, bouncing and the 
transition between them, mass transfer combined with erosion, 
and destructive fragmentation. These are shown schematically in 
Fig. 2, and are discussed in detail in Sects. 3. 1-3.3. In Table 1, we 
give a summary of all the symbols that are used in this section. 



3.1. Sticking and bouncing thresholds 

We consider two dust grains colliding with a relative velocity 
Av. The projectile has a mass m p and the target a mass m, > 
m p . Weidling et al. (2012) found mass-dependent sticking and 
bouncing threshold velocities that can be written as 



Av s tick = I — 



and 



-5/18 



-5/18 



[cms ] 



AVbounce = | — [™S '] 



(3) 



(4) 



where m s = 3.0 ■ 10 12 g and m\, = 3.3 • 10 3 g are two nor- 
malizing constants calibrated by laboratory experiments, and the 



Av 
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proportionality is consistent with the theoretical 



models of Thornton & Ning (1998). The above two thresholds 
mean that collisions with Av < Av st j c k result in 100% sticking, 
and Av > Avb OU nce result in 100% bouncing (provided that nei- 
ther of the particles involved are fragmented). Between these two 
thresholds, we have a region where both outcomes are possible, 
described in more detail in Sect. 3.5. 



Table 1 : Symbols used in the collision model 



Symbol Meaning 



Wp/t 

Av 

Av s tick 

nig 

Avb ounce 

1% 

Vp/t 

M 

/TZrem 

nhm 
Am, 



m 



radius of the projectile/target 

mass of the projectile/target 

relative velocity between the particles 

sticking threshold velocity 

sticking threshold normalizing constant 

bouncing threshold velocity 

bouncing threshold normalizing constant 

center-of-mass velocity of the projectile/target 

relative mass of the largest remnant 

mass of the largest remnant 

mass transferred from the projectile to the target 

accretion efficiency during mass transfer 

mass eroded from the target due to cratering 

net mass change from mass transfer and erosion 

net accretion efficiency from mass transfer and erosion 

velocity needed to fragment with largest remnant fi 

fragmentation threshold normalizing constant 

mass of a monomer (= 3.5 ■ 10~ 12 g) 

total mass of the fragments 



3.2. An energy division scheme for fragmentation 

From the fragmentation with mass transfer experiments de- 
scribed in the previous section, we make the assumption that 
mass transfer at a varying efficiency occurs in all cases where 
the projectile fragments. If also the target fragments, the mass 
transfer is negligible compared to the huge mass loss, and we 
can safely ignore it. We therefore need to determine for each 
collision whether one, both or neither of the particles fragment. 

The majority of the dust collision experiments have however 
been performed between either equal-sized or very different- 
sized particles. In order to interpolate between these two ex- 
tremes, we need to look at the collision energy of the event, and 
determine how this energy is distributed between the two par- 
ticles. Not only the collision energy of an event matters when 
determining the degree of fragmentation, but also the mass-ratio 
between the two particles. In two collisions with equal collision 
energy but with different mass-ratios, we expect the higher mass- 
ratio collision to be less efficient in completely disrupting the 
target, as the energy will be more locally distributed around the 
contact point. In order to take this into account, we choose to 
look at the particles in the center-of-mass frame. In this frame, 
the massive particle is moving more slowly than the small, and 
during the collision moment, the kinetic energy of the particles 
will be reduced to zero. Physically, this corresponds to a fully 
plastic collision where all the energy is consumed by deforma- 
tion and fragmentation. 

We do in this approach assume that the kinetic energy of 
each particle in the center-of-mass frame will be used to try to 
fragment itself. The velocities for the two particles in the center- 
of-mass frame are given by 



Av 



Vn = 



p 1 + m p /m t 
Av 
1 + m t /m p 



v t = 



(5) 
(6) 



All velocities in the center-of-mass frame will from now on be 
denoted as v and will then mean either v p or v t . The above equa- 
tions mean that the largest particle will have the lowest velocity 
in the center-of-mass frame. In the case with an extreme mass 
ratio, m p /m t — > 0, the center-of-mass velocity of the projectile 
and target is given by v p = Av and v t = 0, respectively. 

During a fragmenting collision, the relative size of the largest 
remnant can be described by 



Ttp/t 



(7) 



where m rem is the mass of the largest remnant and ra p / t is the 
original particle mass. Depending on their sizes and material 
strengths, the two original particles can be fragmented to dif- 
ferent degrees. In this model, each collision partner is treated 
individually with a /u t and /i p for the remnant of the target and 
the projectile, respectively. We define the center-of-mass veloc- 
ity required for the largest remnant to have a relative mass \x as 

V 

Blum & Munch (1993) and Lammel (2008) studied the 
threshold velocities needed for two mm-sized particles to frag- 
ment with largest remnants of relative masses ju = 1.0 and 
fi = 0.5, where the former corresponds to the onset of fragmenta- 
tion and the latter corresponds to a largest remnant equal to half 
of the original particle. Beitz et al. (2011) studied the threshold 
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Erosion (E) 

collision velocity 



Fig. 2: Sketch of the five possible outcomes described in Sect. 3 sorted in rough order of collision velocity. Mass transfer and 
erosion act simultaneously in a collision, and we define a mass transfer collision as leading to net growth for the target, and an 
erosive collision leading to net mass loss. This outcome is extremely dependent on the mass-ratio between the particles, thus adding 
a second, vertical dimension to the sketch. 



velocities for cm-sized particles. Interpolating between the re- 
sults for the two sizes, the center-of-mass frame threshold veloc- 
ity can be written as 

v M = (m/m jU )" r [cms 4 ] (8) 

where is a normalizing constant calibrated by the laboratory 
experiments and y = 0.16. The fragmentation threshold velocity 
is given by vi.q, with wi.o = 3.67 • 10 7 g. The velocity required 
for the largest fragment to have half the size of the original par- 
ticle is V0.5, with mo.5 = 9.49 • 10 11 g. The relative mass of the 
largest fragment is fitted to a power law dependent on velocity 
and mass: a 

H(m,v) = C-[^-) ■( T - J ^ T f (9) 
\ 1 g / \ 1 cm s 1 / 

The above equation is valid for all velocities v > vi.o- By fitting 
the fi(m, v) plane to the two parallel threshold velocities given by 
Eq. (8), we get 

a = log(2)/ log(m I . () /mo. 5 ) = -0.068 (10) 
(3 = a /y = -0.43 (11) 
C = m{% = 3.27 [g- a ] (12) 

This means that at a larger collision velocity, the particle will 
fragment more heavily and the size of the largest fragment will 
decrease. More mass is therefore put into the lower part of the 
mass spectrum. 

We can from Eq. 9 determine the largest fragment for each of 
the particles in the collision, and also use it to identify fragment- 
ing collisions. If fi p < 1 and/v, > 1, only the projectile fragments 
and mass transfer occurs. If both fi p < 1 and fi t < 1, both par- 
ticles fragment globally. Since the center-of-mass velocity v is 
inversely proportional to the mass of the particle, we will never 
have a case where only the target fragments and the projectile is 
left intact, even if v p decreases with mass. 



3.3. A new mass transfer and watering model 

We use a new realistic approach to distinguish between colli- 
sions where the target experiences net mass gain due to mass 
transfer, and where it experiences net mass loss due to cratering. 
During each collision, we assume that there is simultaneously: 

- mass added to the target from the projectile via mass transfer 

- mass eroded from the target due to cratering 

We also assume that these two effects act independent of each 
other. This is illustrated in Fig. 3, and can be seen in the 
high-velocity experiments by Teiser & Wurm (2009b). The mass 
change of the largest particle can be described by: 

Am t =m mt -m eT (13) 

where m mt = e ac • m p is the mass added due to mass transfer 
with the accretion efficiency < e ac < 1 and m er is the mass 
lost due to cratering. An increasing velocity not only leads to 
increased mass transfer, but also increased cratering. This makes 
it possible to naturally determine where growth transitions into 
erosion. 

The mass transfer efficiency is obtained from Beitz et al. 
(2011), and depends on both the particle porosity and velocity. 
Since we are unable to track the porosity of the particles, we as- 
sume a constant porosity difference of A<f> = 0.1 between the two 
dust aggregates, where the target always is the more compact 
one. This is likely a reasonable approximation for larger parti- 
cles which have left the hit-and-stick phase and have had time to 
compact during bouncing collisions, which is the region where 
mass transfer can be expected. In our prescription, we also in- 
clude a fragmentation threshold velocity dependence, so that the 
efficiency is always the same for the same degree of projectile 
fragmentation. This results in 

e ac = -6.8 • 10~ 3 + 2.8 • 10~ 4 • Xl£iflE . V (14) 

Vi.o 1 cm s 1 
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m t 



mp-mmt 



trie 




mt-rrier 



Fig. 3: Sketch of the combined cratering and mass transfer pro- 
cess, which occurs during a high-enough velocity collision be- 
tween a projectile with mass m p and a target with mass m t . 
During the collision, an amount m m is added from the projec- 
tile to the target, and the rest of the projectile is put into small 
fragments with a total mass /« p ,f ra g- An amount m a - is simultane- 
ously eroded from the target, and the final mass of the target is 
given by m[ — m t + m ml — m er . 



where vi.o.beitz = 13 cm s _1 is the onset of the fragmentation 
for the 4.1 g particles used by Beitz et al. (2011), and vi.o is 
the fragmentation threshold calculated for the mass of the pro- 
jectile, both calculated using Eq. 8. We here take a maximum 
mass transfer efficiency of e ac = 0.5 as indicated by Wurm et al. 
(2005). Due to the process of fragmentation and mass transfer 
considered here, a higher value would not be reasonable as it 
would point towards complete sticking, which was never ob- 
served at these velocities. 

If the collision energy is not high enough to fragment the 
particle globally, some of the energy is still used to break up local 
bonds between monomers around the contact point, resulting in 
cratering. The cratering efficiency has however only been studied 
in a couple of laboratory experiments. Schrapler & Blum (201 1) 
found for monomer projectiles an erosion efficiency given by 



^ = 1.55-10 ' A! 



m„ 



1 cm s 



-0.4 



(15) 



where m er is the amount of eroded mass and m p = mo is the 
projectile mass, and mo = 3.5 • 10 -12 g is the monomer mass. 
Paraskov et al. (2007) studied the erosion of porous targets both 
with solid and porous projectiles, and found results that varied 
highly with the porosity of the projectile and target. Their results 
are therefore highly uncertain, but roughly agree with an erosion 
efficiency of 

«- 3 Av 

t ^ n I! E I ill e "" ' 



400 1 cm s- 



It should however be noted that for the more compact dust ag- 
gregates that are expected after the compression by bouncing 
phase, the erosion efficiency should be greatly decreased, as is 
generally seen by Teiser & Wurm (2009b). To interpolate be- 
tween the two experiments where the degree of erosion has been 




30 40 
velocity [m/s] 

Fig. 4: The threshold between growth and erosion from the 
model compared to the mass-transfer experiments performed by 
Teiser & Wurm (2009b). Filled circles show experiments where 
the target gained mass, and open circles show where it lost mass. 
The white dotted line shows the threshold from the highly uncer- 
tain erosion prescription from cratering experiments (collisions 
above the line result in erosion, and collisions below result in 
growth). The contours are in intervals of 4% net accretion effi- 
ciency and mark the region with net growth in the final prescrip- 
tion calibrated using the Teiser & Wurm data. 



measured, we assume a mass power-law dependency which can 
be written as 



Av 



1 cm s 



(17) 



where a, b and k are fitting parameters. The above two erosion 
experiments give the efficiency of two different physical effects. 
In monomer impacts, the projectile hits single surface monomers 
and sometimes manages to break the bonds between a couple of 
them. For larger projectiles, restructuring of the target absorbs a 
lot of the collision energy, and a crater is formed both because 
of surface compaction and breaking of monomer bonds. Direct 
comparisons and interpolations between the efficiencies of the 
two effects can not be done without huge uncertainties. A direct 
interpolation between the two effects yields a - 1.55 ■ 10~ 4 , b - 
-0.4 and k = 0.14, but we will below present another way to 
obtain a reasonable erosion prescription. 

As was discussed in the previous section, during a collision, 
erosion and mass transfer usually occurs simultaneously, and the 
net mass change of the target is given by Eq. 13.ForAm t > 0, the 
target experiences net growth, and for Am t < 0, the target expe- 
riences net erosion. With this prescription, the transition region 
is extremely sensitive to the efficiency of the erosion. 

In Fig. 4, we plot the results from the mass transfer experi- 
ments performed by Teiser & Wurm (2009b). We compare this 
to the threshold between growth and erosion (Am t = 0) obtained 
from the mass transfer prescription of Eq. 14 and the erosion 
prescription of Eq. 17. The resulting threshold using the experi- 
ment erosion interpolation is given by the white dashed line, and 
is very pessimistic compared to the mass-transfer experiments. 

Since the experimental erosion prescription is obtained from 
a very different parameter space than what we are interested in, it 
is highly uncertain, and much more so than the mass transfer ex- 
periments discussed below. We therefore choose to calibrate the 
three parameters a, b and k of Eq. 17 using the experimentally 
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Fig. 5: The fragment mass distribution for a 1 cm-sized particle 
after destructive fragmentation events of varying degrees. The 
largest remnant is equal to in a) p = 1, b) p — 0.9, c) p = 0.5, d) 
p = in units of the original particle mass. 



obtained threshold between growth of erosion of Teiser & Wurm 
(2009b). This results in 



1.1 ■ 10" 

-0.4 

0.15 



(18) 



Comparing the net growth efficiency from this fit marked by the 
contours in Fig. 4 to the mass transfer experiments of Kothe et al. 
(2010) (with 1 mm projectiles at velocities of 1-6 m s -1 ) and 
Wurmetal. (2005) (with 1-10 mm projectiles up to 25 m s _1 ) 
results in a rough agreement, even though our model ends up 
slightly pessimistic compared to their results, with net efficien- 
cies roughly half of theirs. Regardless of this discrepancy, we 
take this conservative estimate of the experiments and use it for 
for our model. 



3.4. Fragmentation distribution 

During cratering, mass transfer and destructive fragmentation 
events, the mass of each fragmented particle is divided into 
two parts; the power-law distribution and the largest frag- 
ment. The fragment power-law was determined experimen- 
tally by Blum & Munch (1993) and also used in the model of 
Guttler et al. (2010), and is written 



n(m)dm oc m " dm 



(19) 



where n(m)dm is the number density of fragments in the mass 
interval [m, m + dm], and k = 9/8. 

If the mass of the largest remnant is given by p ■ m, where p 
is the relative size of the largest remnant described by Eq. 9, the 
total mass that is put into the power-law distribution is equal to 
(1 —p) • m. We give the upper limit of the fragmentation distribu- 
tion by min[(l — p),p]-m. This means that as long as p < 0.5, we 
have a single distribution up to the largest remnant. For p > 0.5, 
on the other hand, more than half of the mass is put into the 
largest remnant, which is then detached from the power-law dis- 
tribution. 

This fragmentation recipe is similar to the four-populations 
model of Geretshauser et al. (201 1), with the difference that we 
treat the fragmentation of both particles individually. It is able 
to describe all different degrees of fragmentation, and in Fig. 5, 



we illustrate the fragment distribution for four different values 
of p. In a), we are at the onset of the fragmentation, and all of 
the mass is put back into the remnant, leading to no erosion. In 
b), more than half of the mass is put into the remnant, which is 
therefore detached from the size distrubution, and in c), the ero- 
sion is so strong that that remnant has attached to the power-law 
distribution. Finally, in d), the particle is completely pulverized, 
and all of the mass is put into monomers. 

3.5. Implementation into the code 

In this section, we summarize the conditions and outcome for 
each individual collision type as they have been implemented 
into the code. The different types are, in order, sticking, 
transition from sticking to bouncing, bouncing, mass transfer 
combined with erosion, and destructive fragmentation, and 
they are all shown schematically in Fig. 2. The conditions for 
sticking and bouncing are given in Eqs. 3 and 4, and we use 
Eq. 9 to determine which, if any, of the particles get fragmented 
during a collision, resulting in fragmentation with mass transfer 
or destructive fragmentation. 

Sticking: (Av < Av stick ) 

The two particles stick together and form a bigger particle with 
size Wbig = m t + m v . 

Sticking/bouncing transition: (Av < Av boU nce) 
Transition from 100% sticking to 100% bouncing. We 
assume a logarithmic probability distribution between 
Av stick < Av < Av bounce given by 



1 - ki ■ log 10 (Av) - k 2 



(20) 



where p c is the coagulation probability. We know that at the 
sticking threshold (Eq. 3), the sticking probability is p c = 1, and 
at the bouncing threshold (Eq. 4), the coagulation probability is 
p c = 0. The constants are then 



18/5 



1 log 10 (7w b /m s ) 

log 10 (m p /m s ) 

k-2 = 

log I0 (m b /m s ) 



= 0.40 



(21) 
(22) 



Bouncing: (Av > Av boimce ), (p p > 1) and (p, > 1), (Av < v er ) 
If the collision energy is too high to result in a sticking collision 
but too low to fragment or erode any of the particles, the 
collision results in a growth-neutral bouncing event. The two 
masses involved in the collision are left unchanged. This type of 
collision results in compaction of both particles, but we ignore 
any porosity changes in this model. 

Mass transfer/Erosion: (p p < 1) and (p, > 1) or (m er > 0) 
If the collision velocity is high enough, erosion of the target will 
occur (Eq. 17). Simultaneously, if only the projectile fragments, 
we have a fragmentation with mass transfer event (Eq. 14). The 
resulting mass change of the target is given by Eq. 13. 

The fragmented mass from the projectile is divided into two 
parts, a power-law and the largest remnant, with a total mass of 
m = (1 - e ac )m p . The power-law distribution has a total mass 
of m&ag = (1 - e ac )(l - p P )m P and the largest fragment a mass 
Wrem = (1 ~ fac)y"p'Wp- The fragments excavated from the target 
due to the cratering are also distributed after the power-law 
described in Sect. 3.4, with an upper limit equal to m er . 
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Fragmentation: (jx p < 1) and (p, < 1) 

Finally, if the collision velocity is high enough and the mass- 
ratio is not too large, we get a destructive fragmentation event 
where both particles are fragmented. We treat the fragmentation 
of each particle individually, and get two separate fragment 
distributions, one for the projectile and one for the target. Each 
distribution is divided into two parts; the fragmentation power- 
law distribution with a total mass of mi mg = (1 — yu p / t )m p /t and the 
largest fragmentation remnant with a mass of m rem = /U p /t'« P /t- 

4. The dust-size evolution model 

With the collision model described in the previous section, it is 
possible to describe the collision outcome between any two par- 
ticles. In order to study the evolution of the dust in the protoplan- 
etary disk, however, we need to know about the properties of the 
gas and the sources of relative velocity between the particles. 
In this section, we describe the disk model used in this paper, 
along with the dust evolution code of Birnstiel et al. (2010) that 
has been used together with the new collision model to study the 
dust-size evolution. A summary of the parameters used for the 
disk model is given in Table 2. 

4. 1 . The disk model 

In this work, we follow the dust-size evolution locally at a dis- 
tance of 3 AU from the star. To describe the gas distribution 
over the disk, we use the minimum mass solar nebula (MMSN) 
(WeidenschilHng 1977b; Hayashi et al. 1985). This model is 
based on the current solar system, where the mass of all the plan- 
ets have been used to predict the minimum total mass that would 
have been needed to form them. It however excludes the effect 
of planetary migration and radial drift of dust grains, and the real 
initial disk profile might have been much different (Desch 2007). 
It is however useful for comparison with previous collision mod- 
els. The gas surface density profile of the MMSN is given by 



Table 2: Disk model parameters used in the simulations. 



S g (r) = 1700 



(l Au) 



[g cm 2 ] 



(23) 



where r is the distance to the central star. At 3 AU, this results in 
a gas surface density of 330 g cm -2 , and if we assume an initial 
dust-to-gas ratio of 0.01, a dust surface density of 3.3 g cirT 2 . 

We assume four different sources for the relative veloci- 
ties between dust grains: Brownian motion, turbulence and az- 
imuthal and radial drift. Since we use a local simulation at a set 
point, we take into account the relative velocities that arise, but 
do not allow the particles to move around in the disk. The differ- 
ent sources are discussed briefly below (see Birnstiel et al. 2010 
for a more complete description). 

Brownian motion arises from the thermal movement of the 
particles, and is most effective for the smallest particles. It is 
dependent on the mass of the particles as follows: 



Av BM = 



1 



%k\,T(m t + nip) 



n ■ m t m v 



(24) 



where k\, is Boltzmann's constant, and T = 115 K is the gas 
temperature we assume at 3 AU. 

Turbulent motion arises from the particles interaction with 
the surrounding gas, as it is accelerated by turbulent eddies of 
different size scales. We use the closed form expressions as de- 
rived by Ormel & Cuzzi (2007). The turbulence strength is given 
by the a parameter, generally assumed to lie between 10~ 2 and 



P&rQmctcr 


Symbol 


Value 


TTnit 
UIlll 


Histanpp to star 




3 


AU 


gas surface density 


2 g 


330 


g cirT 2 


dust surface density 


£d 


3.3 


g cirT 2 


gas temperature 


T 


115 


K 


turbulence parameter 


a 


io- 3 




maximum drift velocity 


v„ 


3.9 • 10 3 


cm s" 1 


sound speed 


c s 


6.4 • 10 4 


cm s~ l 


solid density of dust grains 




1.6 


g cnT 3 



10 -4 . The degree at which different particles are affected by the 
turbulence is given by the Stokes number, denoting how strongly 
a particle is coupled to the surrounding gas, which for small par- 
ticles can be written as 



S g 2 



for a < -A mf p 



(25) 



where £ = 1 .6 g cirT 3 is the solid density of the dust grains and 
Eg the surface density of the gas, and A m f p the mean free path of 
the gas. 

Radial drift gives rise to a relative velocity between particles 
as they are differently coupled to the surrounding gas (Whipple 
1972; WeidenschilHng 1977a). This can be written as 

Avrd = |v r (M t ) - v r (M p )| (26) 

where the radial velocity of a particle is given by 

y g 2v n 

Vi = r r (27) 

l+st 2 st + sr 1 

The first term comes from the drag of the surrounding gas on 
the particle as the gas migrates radially, and v g is the the veloc- 
ity of the surrounding gas (Lynden-Bell & Pringle 1974). The 
second term corresponds to the drift of the particle with respect 
to the gas. Due to the gas pressure, the gas moves slightly sub- 
keplerian, and the particle thus experiences a constant headwind 
which causes it to lose angular momentum and drift inwards. 
This effect is strongly related to the coupling between the parti- 
cle and the gas. v n represents the maximum drift velocity, and is 
derived by WeidenschilHng (1977a) as 



Vn = - 



2p g Q k 



= 3.9-10 J cms" 



(28) 



0P„ 



where the gas pressure gradient, p g the gas density and Qk 
the Kepler frequency. 

Azimuthal relative velocities work in a similar fashion to ra- 
dial drift, and arise from gas drag in the azimuthal direction. The 
relative azimuthal velocity can be written as 



Av^ = 



1 



1 



1 +Stf 



i + st 2 , 



(29) 



In Fig. 6, we plot the resulting relative velocity field between 
each particle pair taking into account the four sources described 
above. For particles smaller than ~ 10 yum, Brownian motion is 
the dominant source of relative velocity, causing velocities on 
the order of mrrT 1 . At larger sizes, turbulence becomes impor- 
tant, and velocities quickly increase to ~ lms^'.As can be seen 
in Eq. 27, the radial drift is largest for particles with a Stokes 
number of 1, which at 3 AU corresponds to around 30 cm in 
size. At roughly this size, due to the combined effect of radial 
and azimuthal drift, the particles collide with the smaller parti- 
cles at velocities around 50 m s _1 , which then decreases to 40 m 
s _1 as the particles grow larger and the radial drift decreases. 
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grain size [cm] 



Fig. 6: The relative velocities for each particle pair calculated 
from the four sources described in Sect. 4. 1 using the parameters 
given in Table 2. 

4.2. The Smoluchowski equation 

The dust grain number density n(m, r, z) is a function of the grain 
mass, m, the distance to the star, r, and the height above the mid- 
plane, z, and describes the number of particles per unit volume 
per unit mass. The total dust density can therefore at a point (r, z) 
be written as 

p(r, z) - I n(m, r,z) ■ m dm (30) 
Jo 

and the change in number density with respect to time can then 
be given by the Smoluchowski equation as 

— n(m, r, z) — f f M(m,m',m",r,z) (31) 
ot J Jo 

X n{m! ', r, z)n(m" , r, z) dm' dm" 

where M(m, m', m", r, z) is called the kernel, and describes how 
the mass m is distributed after an interaction between particles m' 
and m". This distribution is determined by the use of a collision 
model like the one developed in this paper and is described in 
Sect. 3. See Birnstiel et al. (2010) on how one constructs M out 
of a collision model. 

In the code implementation, the density distribution is dis- 
cretized over logarithmically spaced mass bins. The resulting 
mass(es) of a collision between two particles will generally not 
coincide with one specific mass bin. In order to solve this, the 
resulting mass is therefore divided between the two neighbour- 
ing mass bins by using the Podolak algorithm described in detail 
by Brauer et al. (2008). 

In order to solve the above equation and track the size- 
evolution of the dust grains, we use an implicit scheme de- 
veloped by Brauer et al. (2008) and Birnstiel et al. (2010). This 
scheme allows for longer timesteps and therefore shorter execu- 
tion times. 

5. Results 

We have performed local simulations of the dust-size evolution 
using the collision model described in Sect. 3 and the evolution 
code briefly described in Sect. 4. In this section, we discuss the 




velocity [cm/s] velocity [cm/s] 



Fig. 7: Comparison between the new collision model (top) and 
the model of Guttler et al. (2010) (bottom). The left and right 
panels show the outcome for equal- and different-sized col- 
lisions, respectively. Green regions mark collisions which are 
growth-positive for the target, yellow marks growth-neutral and 
red marks growth-negative. 'S' marks sticking, 'SB' the sticking 
to bouncing transition, 'B' bouncing, 'MT' net mass transfer, 'E' 
net erosion and 'F', fragmentation. In the transition region, the 
green parallel lines each mark a decrease in sticking probability 
by 20%. 



outcome of the new collision model and compare it to previous 
models. We also show the results of the simulations and compare 
the growth of the large particles to a simple analytical model. 

5.1. The collision outcome space 

With the new collision model, we can determine the outcome 
after a collision between two particles of arbitrary masses and 
velocities. In the upper panels of Fig. 7, we plot the collision 
outcome as a function of projectile size and collision velocity 
for two different mass-ratios. This can be compared to the out- 
come of Guttler et al. (2010) for compact particles shown in the 
bottom panels. It can here be noted that our model naturally de- 
scribes the transition between the two extreme cases of equal- 
sized and different-sized collisions, while Guttler et al. had to de- 
fine a critical mass ratio to distinguish between the two regimes. 
The upper right panel in the figure thus only gives a single snap- 
shot in this transition. 

In the left panels, the two particles are of equal size, and the 
models produce very comparable results. In the new model, the 
sticking region has been enlarged with the inclusion of a transi- 
tion region where both sticking and bouncing is possible. In the 
fragmentation region, the mass-dependent fragmentation thresh- 
old has decreased the velocity needed to fragment larger parti- 
cles, and increased the velocity needed for the smallest particles. 
The end outcome is that the width of the bouncing region has 
decreased significantly. 

In the right panels, the target has a mass 1000 times the 
mass of the projectile. Some important differences can here be 
seen in the fragmentation regime. We can first of all notice the 
new natural transition from growth to erosion which is obtained 
from the balance between growth from fragmentation with mass 
transfer and erosion from cratering (Eq. 13). At this mass-ratio, 
erosion transitions quickly into complete fragmentation. If the 
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grain size [cm] 

Fig. 8: The collision outcome for all pairs of particles with the 
relative velocity field calculated in Fig. 6 and with the same la- 
bels and color code as in Fig. 7. Also included is the net mass 
transfer efficiency, given in intervals of 4%. 



mass-ratio was to be increased further, the fragmentation region 
decreases and is replaced with erosion. 

As long as the projectile is fragmenting, velocities below the 
erosion threshold always lead to growth, and a cm-sized projec- 
tile can initiate mass-transfer already at about 20 cm s _1 (which 
is exactly the result of Beitz et al. 201 1). At Av = 10 m s _1 , pro- 
jectiles smaller than around 1 cm are required for growth. The 
maximum projectile size decreases with velocity as the erosion 
grows stronger, and at Av = 50 m s _1 , growth is only possible 
with projectiles smaller than 100 fim. 

Overall, we predict more fragmentation and cratering than in 
the previous model of Guttler et al. (2010). However, one very 
important change is that growth via fragmentation with mass 
transfer is now possible at higher velocities than the 20 m s _1 
that was the previously predicted threshold, and provided that 
the projectile is small enough, even a collision at 50 m s _1 that 
was predicted in the disk model can lead to growth of the target 
(which was a direct conclusion of Teiser & Wurm 2009b). 

Sticking collisions are also possible at larger sizes, and 
growth-positive mass transfer works at much lower velocities 
than the previously assumed 1 m s -1 . Even if the bouncing re- 
gion shrinks in size, we will see that it will not be enough to 
remove the bouncing barrier. If we insert a particle above the 
bouncing barrier, however, the required relative velocity for it to 
interact beneficially with the particles below the bouncing bar- 
rier has been decreased. These two updates turn out to be quite 
important, as is discussed in more detail in Sect. 5.2. 

The collision outcome for the new model depends on the 
mass of both the projectile and the target, and in the current disk 
model, we use only the average relative velocity between each 
particle pair. This means that a collision between a given pair 
will in the evolution model always result in the same outcome, 
and it would therefore be instructive for us to plot the outcome 
in the particle size-size space. In Fig. 8, we have used the rela- 
tive velocity field calculated in Sect. 4.1 at a distance of 3 AU to 
determine the outcome for each collision pair. 




erosion 



growth 
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Fig. 9: A zoomed-in sketch of the collision outcome space shown 
in Fig. 8. The dashed horizontal line shows the interaction path 
that the seed will experience during its growth. The h parame- 
ter illustrates the minimum distance between the interaction path 
and the erosive region. A positive /; means that the boulder/small 
particle interactions will always be growth-positive, and a neg- 
ative h means that the growth will at some point be stopped by 
erosion. 



In this figure, the bouncing barrier is clearly visible. Dust 
grains of sizes 100 - 800 fim that interact with smaller particles 
will because of too large collision velocities bounce if the parti- 
cle is not smaller than 10 yum. In this case, a small number of col- 
lisions will lead to sticking, but in order to pass the wide bounc- 
ing region, a grain would need to experience 10 9 such sticking 
collisions. The small particles will however be able to coagu- 
late to 100 /im themselves, making growth through the bouncing 
barrier very difficult. 

Collisions between two equal-sized particles larger than 1 
mm will result in destructive fragmentation, but depending on 
what it collides with, a 1 mm-sized particle can also be in- 
volved in sticking, bouncing, mass transfer and erosive colli- 
sions. Because of the fragmentation with mass transfer effect, 
a meter-sized boulder can grow in collisions if its collision part- 
ner is of the right size, in this case smaller than 200 //m. As we 
can see in this plot, the key to growing large bodies is therefore 
to sweep up smaller particles faster than they get eroded or frag- 
mented by similar-sized collisions. 

From Fig. 8, we can already without performing any simula- 
tions see that a cm-sized particle would be capable of growing to 
large sizes if it collides with the right projectiles. The important 
parameter needed to determine this is illustrated in Fig. 9, show- 
ing a sketch of a part of the collision outcome plot. Because of 
the bouncing region, most of the particles will be found in the 
region marked in the figure. A boulder needs to interact ben- 
eficially with these bouncing particles in order to grow, so the 
horizontal interaction path needs to be at all times in the growth- 
positive mass-transfer region. This can be illustrated with the h 
parameter, which gives the minimum difference between the in- 
teraction path and the erosive region. If h is positive, the boulder 
will always interact beneficially with the bouncing particles, but 
if h for some reason would become negative, the growth of the 
boulder would stop. 

We can now point out the interesting effect turbulence has on 
the collision outcome. For particles of sizes between 10 yum and 
10 cm, turbulence is the dominant velocity source. If the relative 
velocity is increased in this regime, the bouncing barrier will 
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Fig. 10: The surface density evolution of the dust population for 
three different simulations at a local simulation at 3 AU. The 
grey diagonal lines correspond to the required surface densities 
for a total number of particles of 1, 10 3 and 10 6 in an annulus 
of thickness 0. 1 AU. In the upper panel, all particles initially 
have a size of 10~ 4 cm, and snapshots are taken between 2 and 
10 6 years. In the middle panel, we have run the same simulation, 
but after 10,800 years, a small number of 1 cm-particles have 
artificially been inserted. In the lower panel, the bouncing barrier 
has been replaced with sticking, allowing the particles to freely 
coagulate to larger sizes. 



be pushed down to smaller sizes. The larger particles are how- 
ever not as much affected by a stronger turbulence, as these sizes 
are also affected by radial and azimuthal drift. This means that 
the h parameter will remain constant or possibly even increase 
with a stronger turbulence. A strong turbulence might therefore 
even be beneficial for this mode of growth, as the larger par- 
ticles will now interact with generally smaller particles, which 
we from Fig. 4 know is beneficial for the mass transfer effect. 
Because of this, even if the boulders due to strong turbulence are 
experiencing relative velocities of ~100 m s , they can grow in 
interactions with the small particles at the bouncing barrier, as 
these have correspondingly decreased in size. 

5.2. The dust-size evolution 

We performed simulations using the new collision model to- 
gether with the local version of the Birnstiel et al. (2010) contin- 
uum dust-size evolution code. In Fig. 10, the mass distribution 
of the particle sizes is given at different timesteps for the three 
different experiments discussed in detail below. 



5.2.1 . Growth up to the bouncing barrier 

In the fiducial case presented in the top panel, the simulations 
are initiated with all dust made up of //m-sized monomers. At 
these small sizes, the relative velocity is driven by Brownian 
motion, and as the particles collide with each other, they stick 
together and form larger aggregates. This leads to a rapid coag- 
ulation phase where the aggregates grow to 100 pm in around 
1000 years. At this point, the particles have grown large enough 
to become affected by the turbulence which quickly increases 
the relative velocities. As we predicted in Fig. 8, the bouncing 
region is too wide to be passed over, and the growth halts at the 
bouncing barrier. 

At this stage, the only particles that can grow are the smaller 
ones, and as time proceeds, more and more particles get trapped 
at the bouncing barrier. This causes the number of small parti- 
cles to keep on decreasing, leading to a continuously narrowing 
size-distribution. After 10 5 years, virtually all particles can be 
found at 100 pm, with very steep tails between 60 and 300 pm. 
If nothing else is done, this is how the dust evolution ends. The 
bouncing barrier efficiently prevents any further growth, and all 
particles remain small. 



5.2.2. A seeding experiment 

In order to investigate the potential of the mass transfer effect, 
we make an experiment where a very small number (i.e. 10~ 1S 
of the total mass) of "seeds" are artificially inserted in the form 
of 1 cm-particles. As can be seen in Fig. 8, the interaction be- 
tween 1 cm and 100 /im-particles results in mass transfer and 
growth of the larger particle, and we expect the inserted parti- 
cles to be able to grow. The seeds are inserted at a single time 
t = 10,800 years, when the particles have reached the semi- 
stable state at the bouncing barrier, and the result can be seen in 
the middle panel of Fig. 10. Exactly how the seeds are formed 
will not be discussed in this paper, but given the small number of 
seeds required, stochastic effects, small variations in local disk 
conditions or grain composition and/or properties might suffice 
to produce them. Some other possibilities are briefly discussed 
in Sect. 5.4. 

In order to better understand the complex interaction be- 
tween all the particles in this experiment that now follows, we in- 
troduce the collision frequency plot given in Fig. 11. This shows 
the collision frequency between each particle pair plotted on top 
of the collision outcome map of Fig. 8, making it possible for 
any given time to identify the dominating collision type for a 
given particle size. 

The first two snapshots in the collision frequency plot are 
taken after 2 and 5900 years, and are identical to the fiducial 
case discussed earlier. At the bouncing barrier, we can see some 
interaction between the 200 pm. particles and the smallest parti- 
cles that do lead to growth, but the frequency is much too small 
to have any significant effect. 

After 10,800 years, the 1 cm seeds are inserted, and they 
grow to larger boulders by sweeping up the small particles 
trapped at the bouncing barrier. As the boulders grow, one can af- 
ter ~200,000 years see a tail of particles with intermediate sizes 
appear behind them. These are formed by the rare collisions be- 
tween the large boulders, and from a single event, two large bod- 
ies have been multiplied to a myriad of fragments also capable 
of sweeping up the particles at the bouncing barrier. This ef- 
fect causes the population of boulders to not only grow in total 
mass, but also in number, which causes a steady and significant 
increase in similar-sized fragmentation. 
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Fig. 11: The collision frequency map for the scenario where 1 cm-particles are artificially inserted at t — 10, 800 years. The interac- 
tion frequency is plotted for each particle pair at six different timesteps plotted on top of the collision outcome space of Fig. 8. This 
makes it possible to identify the dominating interaction for each particle size. Note the always high peak of interactions with parti- 
cles stuck below the bouncing barrier at 1 mm in size. As the large particles grow, they also sometimes collide among themselves, 
producing a tail of particles capable of also sweeping up the bouncing particles. This causes an increase in both mass and number 
for the large particles, and a continuous widening of the size-distribution. 



It can also be noted how the vertical distribution of dust 
around the midplane affects this stage of the evolution. Even if 
there is a huge amount of particles trapped at the bouncing bar- 
rier, they are so small that many of them are pushed out from the 
midplane due to turbulent mixing. The boulders are however so 
large that they have decoupled significantly from the gas, and are 
therefore mostly trapped in the midplane. This causes the sweep- 
up rate of the mm-sized particles by the boulders to be distinctly 
lower than without a vertical structure, and also the internal col- 
lisions between the boulders to be relatively more common. 

The smallest particles that are produced by global fragmen- 
tation and erosion mainly experience two different interactions. 
In the early stages, the smallest fragments are generally being 
swept up by the 100 jum-sized particles stuck at the bouncing 
barrier, since these particles dominate completely both in num- 
ber and mass. They have therefore never any time to coagulate 
to larger sizes themselves, but instead aid in the growth of the 
bouncing particles, which are in turn swept up by the boulders. 
At later stages, as the boulders become more numerous, it is also 
possible that the smallest fragments are swept up directly by the 
boulders. If this growth continues even longer, the two effects 
become equally efficient, and even later, the boulders will start 
dominating in the sweep-up. Regardless of what sizes the small 
fragments interact with, in the end, they are still beneficial for 
the growth of the boulders. 

In the end, a number of 10-70 m boulders have managed to 
form, and the amount of total mass in the large particles have 
increased by the huge factor of 10 12 from what was initially 
inserted into the system, even though the total boulder mass is 
still very small compared to the total dust mass. We find that 
the limiting case for the growth at this point is not so much 



erosion or fragmentation as it is the growth timescale (see also 
Johansen et al. 2008). If the simulation is kept running for longer 
than 10 6 years, the boulders can keep on growing and several 
hundred-meter boulders can form. In other places in the disk 
with higher dust densities and relative velocities, larger boulders 
will be able to form on the same time-scale. 

Growth by sweep-up gives an explanation of how the colli- 
sion part of the growth barrier can be circumvented, but we have 
in these simulations disregarded the effect of the orbital decay 
from gas drag. The growth timescales in Fig. 10 exceed by sev- 
eral orders of magnitude the lifetime of meter-sized bodies sub- 
ject to radial migration. In order to survive, the bodies need to ei- 
ther form on a timescale very much smaller than observed in our 
simulations, which we find unlikely, or there need to exist some 
effect which prevents the orbital decay over an extended pe- 
riod of time (Barge & Sommeria 1995; Klahr & Henning 1997; 
Brauer et al. 2007; Pinilla et al. 201 1). 

5.2.3. Removing the bouncing barrier 

In order to illustrate the importance of the bouncing barrier, we 
make an experiment where all the bouncing collisions are re- 
moved and replaced with sticking. There will therefore be noth- 
ing that prevents the coagulation phase from continuing also to 
larger sizes. The result of this simulation is shown in the lower 
panel of Fig. 10. 

This allows the particles to grow unhindered until they reach 
about 1 mm in size. At this point, they will start to fragment be- 
tween themselves, as we can see in Fig. 8. Since most of the dust 
can be found at this size, heavy fragmentation will occur, and a 
cascade of smaller particles are produced. These small particles 
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grow up to larger sizes again, where most of them again frag- 
ment. From the collision outcome plot, we can however see that 
some particles can be lucky, and instead only sweep up smaller 
particles via mass transfer. They therefore avoid fragmentation, 
and keep on growing, and form the distribution tail extending 
from 1 mm and upwards. 

At a size of roughly 1 meter, the only way to grow is by col- 
liding with 100 pm or smaller particles. But without the bounc- 
ing barrier, most of the mass is instead found in 1 mm particles. 
This means that almost all of the interactions become erosive, 
and the growth halts. 

In the above example, we saw that increasing the stickiness 
of particles in the end actually prevented the growth of large 
boulders. We can from this work draw the conclusion that the 
bouncing barrier might even be beneficial for the formation of 
planetesimals. By hindering the growth of the many, it makes it 
possible for the growth of a lucky few. 

5.2.4. The effect of turbulence 

We have so far studied only one single turbulence parameter of 
a = 10~ 3 . This value is however uncertain for nebular models, 
and in order to investigate the dust evolution for different veloc- 
ity fields, we also study the cases of a — 10~ 4 and a = 10~ 2 . The 
latter represents a strong turbulence which is dominating com- 
pletely over the azimuthal and radial drifts and results in relative 
velocities up to 100 m s . In the former case, the turbulence is 
weak, and contributes very little to the velocities of the larger 
particles. In Fig. 12, we plot the resulting size evolutions for 
the three turbulence parameters but with otherwise identical ini- 
tial conditions. The middle panel is here the same as the middle 
panel of Fig. 10 and is included for reference. 

It can first of all be noted that growth to larger sizes is pos- 
sible even in the case of a very high turbulence, as we predicted 
using Fig. 9. This is due to the decrease in size of the bounc- 
ing particles as the turbulence increases. The smaller impacting 
projectiles can therefore cause the boulder to grow also at higher 
velocities, as seen in Fig. 4. 

The growth timescale is affected by the turbulence in sev- 
eral different ways. Firstly, it increases the relative velocities 
between the particles, leading to higher collision frequencies 
and therefore to more rapid growth. Increased impact velocities 
also increases the net mass transfer efficiency, which in Fig. 4 
is seen to be particularly low for velocities below 10 m s _1 . 
An increased turbulence also mixes small particles further out 
from the midplane, which decreases their midplane densities 
where the largest boulders gather, lowering the growth rate. 
From Fig. 12, it is clear that the growth-positive effects are dom- 
inating, leading to enhanced growth with increased turbulence. 
In the low turbulence case, the growth is especially slow for the 
sizes between 1 cm and 100 cm. In this regime, the relative ve- 
locities are very low since the contributions from the azimuthal 
and radial drift are small, causing very low net mass accretion 
efficiencies during the sweep-up growth. 

In the a = 10~ 2 case, we can in the final snapshot see that 
a separate peak has appeared for the largest boulders. It forms 
when the boulders have fragmented so much between them- 
selves that the large intermediate-sized fragments are at num- 
ber densities roughly equal to the boulders, which in this case 
happens after ~ 600, 000 years. At this point, the boulders are 
significantly eroded and fragmented, which creates even more 
intermediate-sized particles capable of yet more fragmentation. 
This results in a fragmentation cascade, and a rapid flattening of 
the size-distribution. This does however not cause all the large 
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Fig. 12: The surface density evolution of the dust population for 
three different turbulence strengths (a = 10~ 4 , 10~ 3 , 10~ 2 ) at 
3 AU, with seeds inserted artificially after t = 10, 800 years. 



particle to be destroyed, and those that survive can continue 
to grow by sweeping up the small bouncing particles that still 
dominate both in number and total mass. This effect occurs also 
for the cases of weaker turbulence if the simulations run long 
enough, but is never severe enough to cause a complete halt of 
the growth. 

We have have here shown that growth can proceed also in re- 
gions of high turbulence. An MRI turbulent disk might however 
also create an additional velocity source, where local gas den- 
sity fluctuations can excite the orbital eccentricities of the plan- 
etesimals. Ida et al. (2008) found that this can lead to velocities 
beyond break-up even for small planetesimals, but only studied 
collisions between equal-sized bodies. In our model where only 
the material strength of the body is included, a collision between 
two large and similar-sized bodies will be destructive even in re- 
gions of very low turbulence, so in this regime, it will affect our 
conclusions very little. Excited planetesimal orbits will however 
also increase the impact velocities of the small particles, which 
could pose a problem. In the previous cases, growth has been 
possible in regions of high turbulence because the sizes of the 
bouncing particles have simultaneously been suppressed by the 
turbulence, but this is not necessarily the case in the scenario 
studied by Ida et al. This might cause an increased erosion of the 
growing planetesimals, but investigating this further is outside 
the scope of this paper. 
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5.3. A simple analytical model 

To verify what we see in the simulations, where growth is caused 
by the sweep-up of small particles by a few lucky seeds, we will 
make a simple analytical growth model in a similar fashion to 
Johansen et al. (2008) and Xie et al. (2010). We consider a sce- 
nario where a single large body of mass m is moving around in 
a sea of small particles of mass m s stuck at the bouncing barrier. 
The growth of the large body can then be described by 



dm 

— = crAv ■ e net n s m s 
at 



(32) 



where cr - k (a s + a) 2 ~ na 2 is the collisional cross-section, Av 
the relative velocity, n s the number density of the small particles 
and e net the average net mass gain efficiency of a collision. The 
change in size a with respect to mass can be written as dm = 
An^a 2 da, where £ = 1.6 g crrT 3 is the internal density of the 
large body. We can now write 



da 1 p s 

= Enet A V 

df 4 £ 



(33) 



where p s = n s m s is the mass density of the small particles. We 
will for now assume that the relative velocity is caused only by 
turbulence, and can therefore for small particles be estimated by 
(Weidenschilling & Cuzzi 1993; Cuzzi et al. 2001) 



Av = ^l-aSt-c s 



(34) 



where a = 10~ 3 is the turbulence parameter, St the Stokes num- 
ber for the boulder given by Eq. 25 and c s the sound speed in the 
gas given by 



c s = 




(35) 



where k\, is Boltzmann's constant, T = 1 15 K is the temperature 
of the gas, p. - 2.3 the mean molecular weight and m p the mass 
of a proton. 

We assume p s to be constant and unaffected by the sweep-up 
of the large boulder. When the boulder has grown large enough, 
it will be very little affected by the vertical mixing from turbu- 
lence, and will mostly be found near the midplane. We therefore 
take p s = 10" 12 g cirT 3 to be the midplane density assuming a 
Gaussian vertical profile with a surface density X s = 3.3 g cm -2 
and that all the of the dust mass can be found in the small parti- 
cles. Eq. 33 can then be solved analytically: 



a(f) = 



1 A 9 € 



pfnety -J^^-TTff ' C s (f- fo) + V^O 



(36) 



In Fig. 13, we plot the above analytical solution to the boulder 
growth assuming e net = 0.006 which is the expected efficiency 
for the interaction between a 1 cm and a 1 mm body. We as- 
sume the other parameters identical to the simulations. This is 
compared to the observed growth of the largest particle from 
the seeding simulations in Sect. 5.2.2 using the full collision 
model. As can be seen, our analytical solution compares well 
to the growth rate of the simulations for the first stage of growth, 
and then at later stages differ greatly. This is mainly due to two 
reasons. Assuming a constant e net means that it at larger sizes 
will be underestimated, as it in the full model depends on the 
collision velocity and varies between and 0.15. The relative 
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Fig. 13: The growth of a boulder in the simulations (black) com- 
pared to two growth estimates. The red line corresponds to the 
simplified analytical estimate, and the blue line to the estimate 
using the full relative velocity and accretion efficiency prescrip- 
tions. 



velocity prescription in Eq. 34 only takes into account the turbu- 
lence, and is also not valid for large particles. In order to take this 
into account, we use the full net growth efficiency from Eq. 13 
and all four velocity sources discussed in Sect. 4. We then solve 
Eq. 33, and plot the resulting solution in Fig. 13. The new so- 
lution and the simulation results now agree very well, only to 
differ slightly at the largest particle sizes. 

We find that this enhanced growth is triggered by the exis- 
tence of the tail of intermediate-sized particles that form from 
fragmenting events between large bodies. The largest boulders 
found in the simulations are those that manage to grow by avoid- 
ing interactions with other large bodies, and if they manage to 
interact beneficially with some of the intermediate-sized frag- 
ments, the growth rate is enhanced. 

5.4. Forming the first seeds 

The width of the bouncing barrier has with recent sticking and 
fragmentation experiments decreased, but is still solid enough 
to prevent particles from growing through it. As previously dis- 
cussed, this might even be positive for the planetesimal forma- 
tion, as it prevents too many large bodies from forming, keeping 
them from fragmenting and eroding between themselves. It is 
however necessary for some cm-sized seeds to be formed that 
can initiate the sweep-up. 

We have in this work found that the bouncing barrier can 
be circumvented, and growth via mass transfer initiated, even if 
only a very small number of cm-sized particles is introduced to 
the system. We have in this paper not investigated in detail how 
this might happen, but there are many possibilities. Larger par- 
ticles may form outside the snow-line mixed with ices, and drift 
inwards to a region where sweep-up becomes possible. Calcium- 
and alumunum-rich inclusions (CAIs) are cm-sized particles that 
are believed to have been formed early near the Sun and trans- 
ported outwards in the disk (Ciesla 2009), and could also make 
up the first generation of seeds. 

It is also possible for some lucky particles to grow through 
the barrier simply by interacting with enough smaller particles, 
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as can be seen in the collision outcome plot in Fig. 8. In or- 
der for this to happen, however, it is necessary to have a wide 
enough size-distribution so that there is enough of the small par- 
ticles that the lucky particles can interact with. This is difficult, it 
turns out, as we have in these simulations found that the particle 
size-distribution will quickly narrow as all the small dust coag- 
ulate up to the bouncing barrier. If extra mass is continuously 
introduced to the system the first 10,000 years or so, for exam- 
ple due to nebular infall of the collapsing protostellar cloud, the 
particle size-distribution could be wide enough for some seeds 
to be formed. Beitz et al. (2012) also found in laboratory experi- 
ments that chondrules and dust aggregates tend to stick at higher 
velocities than collisions between two dust aggregates, so that 
early chondrules could grow where dust could not. 

Another possibility is the introduction of a velocity distri- 
bution for each particle-size. Most work with the dust-size evo- 
lution consider only the average relative velocity between two 
particles, but in reality, some particles will collide also at much 
higher and lower velocities. If some particles are lucky enough 
to only collide with others at low velocities, they might stick 
together where they would otherwise only bounce, and in that 
grow large enough to initiate sweep-up. 

6. Discussion and Conclusions 

In order to explore the possibility of growth to planetesimal 
sizes by dust collisions, we have implemented a new collision 
model which is motivated both physically and experimentally. 
It is streamlined to work with continuum codes, while still able 
to take into account all important collision types that have been 
identified in previous work. Using it with continuum dust evolu- 
tion codes makes it possible to resolve all dust grains regardless 
of their numbers, something which we have found to be essential 
for the study of growth of dust grains above cm-sizes. 

Even though collisions between large similar-sized dust 
grains generally lead to fragmentation, this is not necessarily 
true for larger mass-ratios if the projectile is small enough. As 
shown in laboratory experiments, if the projectile is smaller than 
0.1-1 mm in size, fragmentation with mass transfer can cause 
growth of the target at impact velocities as high as 60 m/s. In 
order for dust grains to grow larger than cm-sizes via collisions, 
the number of large particles therefore has to be very small to 
avoid destructive fragmentation among themselves, and most of 
the dust mass kept at small sizes. 

We have found in our simulations that direct growth of plan- 
etesimals via dust collisions still is a possibility. We find that 
the bouncing barrier introduced by Zsom et al. (2010) might be 
beneficial or even vital for the planetesimal formation, as it pro- 
vides a natural way to keep most of the dust population small. 
This small dust makes ideal material for the sweep-up process 
if larger bodies manage to form. By artificially inserting a few 1 
cm-sized seeds in our simulations, we find that they can sweep 
up the small dust via fragmentation with mass transfer and grow 
to ~ 100 meters in size on a timescale of 1 Myr. This opens up ex- 
citing new possibilities that need to be taken into account when 
studying dust growth in protoplanetary disks. 

The growth rate is in our simulations relatively slow, mainly 
due to a low mass transfer-efficiency, the high turbulence kicking 
the small particles away from the midplane where the boulders 
are concentrated, and the low dust densities already at 3 AU. This 
means that it might be necessary to form planetesimals in regions 
with enhanced densities such as in pressure bumps, where it can 
be accelerated, and also that there is a need for the radial drift 
to be prevented over long timescales. We find in this work that 



planetesimals formed in this way would be small, between 100 
meters to some kilometers in size, which is smaller than what 
was generally predicted by Xie et al. (2010), where a larger mass 
transfer-efficiency was assumed. 

This mode of growth is however able to function also in high 
turbulence regions with relative velocities reaching 100 m s . 
This is because the turbulence simultaneously suppresses the 
size of the bouncing particles, which are therefore more likely 
to cause mass transfer rather than erosion of the boulders during 
impacts. Another important effect we have found is that a very 
small number of seeds are necessary to initiate sweep-up. As the 
large particles collide with each other, they will create a num- 
ber of intermediate-sized particles that are also able to sweep-up 
the small particles, causing the population of large particles to 
increase not only in total mass, but also in number. 

Exactly how these seeds would be introduced into the sys- 
tem needs to be explored further. We have shown that thanks to a 
combination of more efficient sticking found in recent laboratory 
experiments and our ability to numerically resolve very small 
numbers of particles, the bouncing barrier can be overcome or 
circumvented. If the dust size-distribution is wide enough, a 
small number of lucky particles might grow over the barrier via 
hit-and-stick collisions with much smaller particles. This shows 
that the initial conditions of planet formation might be very im- 
portant for how the dust growth proceeds. Another option, which 
we intend to investigate further, is the effect of adding a velocity 
distribution for each particle size, which would make it possi- 
ble for some lucky particles to always experience low-velocity 
collisions and thus grow to large enough sizes. 

This work shows where the focus of future laboratory exper- 
iments should lie. It is clear that in order for collisional growth 
of larger particles to be possible, it has to occur between parti- 
cles of very different sizes. However, very few experiments have 
been performed to quantify the amount of erosion or mass trans- 
fer for varying projectile sizes, porosities, velocities and impact 
angles, and the understanding of the physical process is still not 
clear. The maximum size a projectile can have to cause growth 
of the target due to mass transfer determines whether the small 
bouncing particles cause erosion or growth of the large boulders. 
It is also necessary in the lab to determine the smallest projec- 
tile size that still leads to growth. Laboratory experiments show 
that monomer impacts lead to erosion, but whether erosion still 
occurs at 10 or 100 fxm is not known. In order to determine in 
which parameter space fragmentation with mass transfer occurs, 
more experiments need to be performed. 

With an implicit scheme, this code runs fast in terms of ex- 
ecution time, and should be suitable for global disk simulations, 
something which is not possible with a Monte Carlo approach. 
In future work, we intend to implement the new collision model 
in the global dust evolution code of Birnstiel et al. (2010). This 
will make it possible to naturally study the number of seeds that 
can be formed in the disk and how they migrate inwards to other 
parts of the disk, and also to see how this more sophisticated 
collision model affects the dust-sizes further out in the disk at 
80-100 AU where comparison with observations is possible. 
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